Analysis of quantum coherence in bismuth-doped sihcon: a system of strongly coupled 
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We investigate bismuth-doped silicon (Si:Bi) in the "nuclear-electronic" regime where the hyper- 
fine coupling of 1.4754 GHz is comparable to the electron Zeeman term in the spin Hamiltonian. 
This corresponds to intermediate magnetic fields, B ~ 0.1 — 0.6 T. There is growing interest in Si:Bi 
as an alternative to the well-studied proposals for silicon based quantum information processing 
(QIP) using phosphorus-doped silicon (Si:P). We focus here on the implications of the anomalously- 
strong hyperfine coupling. In particular, we analyse in detail the regime where recent 4 GHz pulsed 
magnetic resonance experiments have demonstrated orders of magnitude speed-up in quantum gates 
involving transitions which are forbidden at high fields. We also present calculations using a phe- 
nomenological Markovian master equation which models the decoherence of the electron spin due 
to Gaussian temporal magnetic field perturbations. The model quantifies the advantages of certain 
"optimal working points" identified as the = regions, which come in the form of frequency 
minima and maxima. We show that at such regions, dephasing due to the interaction of the electron 
spin with a fluctuating magnetic field in the z direction (usually adiabatic) is completely removed. 

PACS numbers: 03.67.Lx,03.67.-a,76.30-v,76.90.-|-d, 
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I. INTRODUCTION 

Beginning with the seminal proposal by Kane there 
has been intense interest, for over a decade, in the use 
of Si:P [2[ as qubits for quantum information processing. 
This donor-impurity spin-system continues to demonstrate 
an ever-increasing list of advantages for manipulation and 
storage of quantum information with currently available 
electron paramagnetic resonance (EPR) and nuclear mag- 
netic resonance (NMR) technology. The Si:P system has 
four levels due to the electron spin S — ^ coupled to a ■^^P 
nuclear spin I = 

The key advantages are the comparatively long decoher- 
ence times, which have been measured to be of order mil- 
liseconds for the electron spin, for natural Si:P. They are of 
order seconds for the nuclear spin, so the nuclear spin has 
been identified [l| as a resource for storing the quantum 
information. For all but the weakest magnetic fields (i.e. 
Bo > 200G) i, the electron and nuclear spins are uncou- 
pled so may be addressed and manipulated independently 
by a combination of microwave (mw) and radio-frequency 
(rf) pulses respectively. The two possible electron-spin 
transitions correspond to EPR spectral lines, while the nu- 
clear spin transitions are NMR lines. Nuclear spin-flips are 
much slower: a 7r-pulse in the NMR case takes three orders 
of magnitude longer than for the EPR-allowed transitions. 

However, over the last year or so, there has also been in- 
creasing interest in another shallow donor impurity in sil- 
icon, the bismuth atom 043- The SiiBi system is unique 
in several respects: it is the deepest donor, with a binding 
energy of about 71 meV; it has a very large nuclear spin, 
I = 9/2; it has an exceptionally large hyperfine coupling 
strength, ^ = 1.4754 GHz; in addition, ^"^Bi is mono- 
isotopic. Recent measurements of the decoherence times in 
natural silicon have revealed T2 times of order 30% larger 
than for Si:P, an effect attributed to the smaller Bohr ra- 
dius of Si:Bi Q. The dominant decoherence process is the 



spin diffusion [l3| associated with the / = 1/2, ^^Si nu- 
clei occupying just under 5% of sites in natural silicon; the 
dominant ^*Si isotope has no nuclear spin and thus does 
not contribute to the dipole-coupled flip-flop process which 
drives the spin diffusion. A recent study of P donors in ^^Si 
purified to such a high degree (less than 50 ppm of ^^Si) 
that spin diffusion may be neglected, revealed T2 times po- 
tentially up to 10s [ll|. Although studies of isotopically en- 
riched Si:Bi have yet to be undertaken, decoherence times 
of the same order are possible. The coupling with ^^Si was 
investigated in The very large nuclear spin I ~ 9/2 and 
associated large Hilbert space may provide a means of stor- 
ing more information Q • Efficient hyperpolarization of the 
system (to about 90%) was demonstrated experimentally 
in i. 

The present study investigates the implications of the 
very large hyperfine coupling, ^ = 1.4754 GHz of Si:Bi, 
as well as its large nuclear spin. While mixing of the Zee- 
man sublevels \ms,mi) is not unexpected and has even 
been investigated for Si:P for weak magnetic fields Q, for 
Si:Bi the regime where the hyperfine coupling competes 
with the external field, i.e. A ^ fiBo/h, is attained for 
magnetic fields Bq ~ 0.1 — 0.6 T which are moderate, but 
within the normal EPR range. In a previous paper [7[ we 
identified interesting consequences in this range of mag- 
netic fields. In particular, wc identified a set of points and 
field values associated with "cancellation resonances" , thus 
named because of an analogy with the electron spin echo 
envelope modulation (ESEEM) phenomenon of "exact can- 
cellation" [T^]. At the cancellation resonances, as we see 
below, the system Hamiltonian takes a simpler form. Here 
we investigate the relative advantages of these regimes for 
QIP. 

A recent experimental study using an S-band (4 GHz) 
pulsed EPR spectrometer [l^ identified a set of 4-states of 
Si:Bi which are, at high fields, entirely analogous to the 4- 
level subspace of Si:P but have the advantage that all four 
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possible transitions may be driven by fast EPR transitions 
(on a nanosecond timescale), in contrast to Si:P, where the 
two NMR transitions are three orders of magnitude slower. 
The study showed that for Si:Bi in the intermediate field 
regime, it is possible to equalize the time for a 7r-rotation 
of two transitions which correspond to an EPR and NMR 
transition respectively, at high fields. We show that this 
represents a novel scheme to implement universal quantum 
computation using, exclusively, fast EPR technology. This 
represents a clear advantage offered by the Si:Bi system. 

In 01 , it was suggested that the "cancellation resonance" 
points might offer further advantages in minimising and 
controlling exchange coupled perturbations. Here we do 
not investigate this, but consider effects of Z-noisc and 
X-noisc (in other words Gaussian temporal magnetic-field 
fluctuations along the X and Z direction) on dccohercnce. 
This may be relevant to the behavior of isotopically en- 
riched Si:Bi. We show that for Z noise, which is usually 
adiabatic, the = points offer decoherence-free zones. 
In analogy with work done on superconducting qubits [l^ , 
wc call these "optimal working points" . The system docs 
not show such advantages for X noise, however, which leads 
to temperature-independent depolarising noise. 

In section |TT] we introduce the Hamiltonian for coupled 
nuclear-electronic spin systems for S = 1/2 and arbitrary 
/. We explain the analogy between the EPR hyperfine can- 
cellation resonances and the Electron Spin Echo Envelope 
Modulation (ESEEM) regime of exact cancellation. We 
also introduce the individual resonances, describe salient 
features and present time-dependent calculations. In Sec- 
tion |lTT] we introduce the system as a pair of coupled qubits 
and compare with Si:P. We propose here a scheme of uni- 
versal two-qubit quantum computation. In section ITVl we 
introduce a model of decoherence caused by a temporal 
fluctuation of the external magnetic field, and study the 
effect of the cancellation resonances on the decoherence 
rates this model predicts. We conclude in Section fVl 

II. THEORY OF COUPLED 
NUCLEAR-ELECTRONIC SPIN SPECTRA 

A. The Hamiltonian 

nuclear-electronic spin systems such as Si:P and Si:Bi are 
described by the Hamiltonian: 

Ho=LOoSz-LJoSiz+AS.i (1) 

where ujq represents the electron Zeeman frequency given 
by Bg/3. Here B is the strength of the external magnetic 
field along the z direction, g is the electron g-factor, and /3 
is the Bohr magneton. 6 ~ uji/ojq = 2.488x 10~^ represents 
the ratio of the nuclear to electronic Zeeman frequencies. 
A is the isotropic hyperfine interaction strength. The op- 
erators S and I act on the electronic and nuclear spins, 
respectively. The total spin of the electron and nucleus is 
given by S and /, and their z-axis projection by 7715 and 
m/, respectively. 

For the systems considered, the electron spin is always 
S = ^. As a result the dimension of the Hilbert space is 
determined by the particular nuclear spin: for a given nu- 
clear spin /, there are 2(2/ -I- 1) eigenstates which can be 
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FIG. 1: The 20 spin energy-levels of Si:Bi may be labelled in 
order of increasing energy |2)...|20). States |10) and |20) are 
not mixed. State |10) is of especial significance since it, rather 
than the ground state, is a favourable state to initialise the sys- 
tem in. Experimental hyperpolarization studies concentrate 
the system in this state. Thus, in our coupled 2-qubit scheme, 
state 1 10) corresponds to our |OeO„) state; in the same scheme, 
states |9) = |Oeln) and = |leOn) are related to |10) by a 
single qubit flip while for |12) = [leln) both qubits are flipped. 

superpositions of spin basis states jm^, mi) . However, since 
[Hq^Sz + Iz] = 0, the Hamiltonian in Eq.(IT]) decouples to 
a direct sum of one and two-dimensional sub-Hamiltonians 
and with constant m — rris + mj. The former act 
on the basis states \nis = ±i,m/ = ±(/ -I- i)), while the 
latter act on the basis states \ms = ±^,mi = m =F i) such 
that |m| < /+ i. The two-dimensional sub-Hamiltonians 
can be expanded in the Pauli basis. In particular, the ex- 
ternal field part of the sub-Hamiltonian operator is given 
in matrix form: 

uoSz - ujoSiz = Y [(1 + 6)cTz - 2mSl] . (2) 
The z-componcnt of the hyperfine coupling: 

ASz®iz^^{maz~t/2) (3) 

is seen to have an isotropic component as well as a non- 
isotropic component dependent on az, while the x, y com- 
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ponents are given by: 

A{S, (E}t+Sy® iy) = ^{1(1 + 1) + ^ - (4) 

Summing the above terms gives each H^, whereas only 
Eqs. © and © contribute to H^^f: 



TT2d _ 



rrld 

-"m=±(/+i) 
^rj7 



= m + CJo(l + 5) 

= (/(/ + !) + i-m^/^ 
= i(l + 4wom5). 



(5) 



Wo = ^ is the rescaled Zeeman frequency. We define a pa- 
rameter i?^ = + where Rm represents the vector 
sum magnitude of spin x and z components in the Hamiho- 
nian. We denote 9m as the inchnation of Rm to the z-axis. 
such that cos 0m = ^m/Rm ^nd sin 0m = ^m/Rm- Thcn, 
H?f can also be written: 



H. 



2d 



A 



{Rm cos 0m<Jz + Rm sin 6m(Tx - Cml) (6) 



Em=±(I+l/2) 



(11) 



It is important to stress that the Cz, ax above are quite 
unrelated to the 5^, electronic spin operators. They are 
simply a method of representing the two-dimensional sub- 
Hamiltonians. It also makes obvious that the magnetic 
fields for which m + ujo{1 + 5) = are special: here the 
term in H^^ dependent on cr^ vanishes entirely. These cor- 
respond to Landau-Zcncr crossings. The point at which 
m + ujo{l + 6) = for 7Ti = — (/ + i) also has special prop- 
erties, although different than those of the previous case. 
In 01 a further point specific for Si:Bi, for which wq = 7 A, 
was also identified; here oc ((Tx + fz) (ignoring the triv- 
ial term proportional to the identity). These field values 
were all classed as "cancellation resonances" and all have 
potential applications in QIP. 

In Figlllthe exact expressions in Eas. (fTO)) . (|TT|) were used 
to reproduce the spin spectra investigated for Si:Bi in e.g. 
[3| and @. These equations can be used to describe any 
arbitrary coupled nuclear-electronic spin system obeying 
the Hamiltonian Eq.([l]), such as other donor systems in Si 
including P and As. However, throughout this paper we 
only use Si:Bi as the instantiated system for which we give 
numerical solutions. 



The range of values that 0m can take are given by: 



B. EPR transitions 



[0,arctan(^|^ 
[0, ^ -I- arctan 



|m| 



such that 771 > 
such that m ~ 
such that 771 < 



(7) 



The EPR transitions |-|-, 77i) -r-> | — , tti — 1) are dipole al- 
lowed at all fields. Their relative intensities are 
/ « 2/+^-_J(m, = ^\Sx\m', = Since - 

i) P = 1/2, variations in line intensities arise from the mix- 
ing of the states. 



where the minimal value occurs as _B — )■ oo and the maxi- 
mal value occurs at B = 0. Note that 9m < tt V B. 

Straightforward diagonalisation of then gives the 
eigenstates at arbitrary magnetic fields: 



where 



cos ( ^ ) , 



1 J. Ix 



±1 " III 
sm — 
2 



(8) 



(9) 



and with the corresponding eigenenergies: 



-i(l+4wom(5)±i?„ 



For H^, 9m ~ 7r/2 whenever m + ojq{1 + 6) ~ which 
occurs only for < m < 0. The B ^ oo limit cor- 

responds to 9m 0, such that the eigenstates in Eq.® 
tend to the \ms,rni) basis states, which we call the high- 
field limit states. The higher the hyperfine constant A, the 
higher the magnetic field would have to be for these eigen- 
states to be approximated by the high-field limit states. 

has 9 m = Oy m, and hence gives the uncoupled eigen- 
states I ± ^,±(/ + 5)) at all magnetic fields. These have 
the simplified eigenenergies: 



1 



(1 + COS0,„)(1 + COS^m-i) 

(12) 

If mixing is significant |+,r77.) •H- |+,?Ti — 1) transitions (of 
intensity /+) and | — ,?n — 1) <-> \ — ,m) transitions (of in- 
tensity /~), which are EPR-forbiddcn (but NMR-allowcd) 
at high field, become strong, with relative intensities: 



-^m->m-l I'^m.l l^m-ll 



and 



-(1 -I- C0S6',„)(1 - COS^m^i) 

(13) 



m_l_,m OC |a„_info„J^ = i(l - C0s6'm)(l +C0s6'm_i) 



(14) 

EPR-forbiddcn lines disappear at high fields; as wq — > 00, 
(10) one can see from Eq.® that /±,m^m-i ^ ^ since 

l^mP oc ^ at high fields. Eqs.(Iini),(IIll), (US]) and (HH) are 

exact so are in complete agreement with numerical diago- 
nalisation of the full Hamiltonian. 



C. The uio — 1,2,3,4 cancellation resonances and the 
frequency minima 

These cancellation resonances are Landau-Zener or 
avoided crossings, when the az term of vanishes, lead- 
ing to "Bell-like" eigenstates in the spin basis. In FigI2]we 
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FIG. 2: (a) Comparison between theorv(Eas. (fTO)l . (fTTjl and 
p2|l ') (black dots) and experimental CW EPR signal (red on- 
line) at 9.7 GHz. Resonances without black dots above them 
are not due to Si:Bi; The large sharp resonance at 0.35 T is 
due to silicon dangling bonds while the remainder are due to 
defects in the sapphire ring used as a dielectric microwave res- 
onator. The variation in relative intensities is mainly due to 
the mixing of states as in Eq.([8|. The variability is not too 
high but the calculated intensities are consistent with exper- 
iment and there is excellent agreement for the line positions, 
(b) Calculated EPR spectra (convolved with a 0.42 mT mea- 
sured linewidth); they are seen to line-up with the experimental 
spectra at / = ^ = 9.7 GHz). The "cancellation resonances" 
are indicated by arrows and integers —m = 0,1,2,3,4.... The 
Ldo = 7 A maximum of ^ indicated by the integer 7 corresponds 
to that shown in the < 2 GHz Electron Nuclear Double Reso- 
nance (ENDOR) spectra of [||. 



see several transitions of type |+,™) ^ | — ,m — 1) that 
have a minimum frequency, closely related to the cancella- 
tion resonances. We can show that these minima (obtained 
by setting ^'-"'^"'g^"''^-' = 0) occur for: 

cos 9m = — cos 6m-l (15) 

if 7Ti < 0. This implies: 



which means that the minima lie exactly midway in angular 
coordinates between cancellation resonances. For example, 
for the |12) o |9) line (| + ,m = -3) o |-,m = -4)) the 
minimum is at wq = ^ = 3.57A so B = 0.188 T. Here, 
i?_4 has passed its resonance point at 0.21 T (for which 
6'_4 = 7r/2) by an angle = arccos and R-z is at an 
equal angular distance before its resonance at ~ 0.16 T. It 
has already been shown that a reduction in df/dB results 
in experimentally measured reductions in decoherence time 
0. The df/dB = points imply a reduction in the sensi- 
tivity to temporal magnetic fluctuations. In section ITVl we 
study how these frequency minima affect decoherence rates 
caused by Gaussian temporal magnetic-field fluctuations. 

D. The LjQ — 7 resonance; a frequency maximum 

The cluster of lines in Figl2jb) at frequencies / < 1 GHz 
represent |±,m) — > |±,m— 1) transitions EPR-forbiddcn 
at high B (their intensity is seen to tend to zero for B > 
0.5 T). These transition frequencies have maximal values 
calculated by setting — which leads to : 

cos6',„ = cos6'„_i (17) 

if TO < 0. Since Q < 0„i < tt, this implies that 
= Om-i- Only one such maxima, for the transition 
|±,-3) o |±,-4) occurring at (Do = 7 and B ~ 0.37 T 
can be observed by EPR, however. The other maxima oc- 
cur at fields i? > 0.5 T for which the EPR line intensities 
become vanishingly small. The wq = 7 frequency minima 
is also classified as a cancelation resonance because of the 
symmetric property that the involved exhibit. This 
resonance occurs when 0_3 = = 7r/4, which implies 
that H'i% oc Hli (x{aa:+<Jz)- 

The ujQ — 7 cancellation resonance offers possibilities for 
more complex manipulations. It has been suggested [1, Q 
that the larger state-space of Si:Bi may be used to store 
more information. Thus we can show that at ujq — 7 A, 
a single EPR 80 ns) pulse can map any coherences be- 
tween the m = —4 states into the same coherences between 
the m = —3 states. The condition Eq. pT)) implies that the 
amplitudes a^g = and 6^3 = ±6^4. This means that 
an EPR pulse will effect the rotations |12) |11) and 
|9) |8) at the same rate. For instance, if the initial 
two-qubit state is I*!*) = cii|12) -I- C9I8) a 7r-pulse will yield 
l^*) = cii|ll) — C9I9), and so produces a mechanism for 
temporarily storing the two-qubit state (within a relative 
TT phase shift). This is illustrated in FiglS) 

E. The ujo = 5 resonance : two-photon transitions 

The TO — —5 state, |10), is not associated with an anti- 
crossing at any field as the Hamiltonian leaves it uncoupled 
to any other basis state. Nevertheless the fields for which 
m + ujQ w (at ujQ ~ 0.26 T) represent the most dras- 
tic type of cancellation: the A_5 term in fflg vanishes. 
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ences between state |10) and both states |11) and |9). In 
effect one may use two-photon, second order processes to 
transfer population between states |9) and |11) (recall that 
simultaneous spin flips are forbidden for isotropic hyperfine 
coupling). Fig|3] illustrates this. 



F. Analogy with "exact cancellation" 

Exact cancellation is a widely used "trick" in ESEEM 
spectroscopy. A coupled nuclear-electronic system with 
anisotropic hyperfine coupling, which is weak compared 
with electron spin frequencies (on the MHz scale rather 
than GHz scale), has a rotating frame Hamiltonian [l^ : 



FIG. 3: Shows that near the lio = 7 cancellation resonance, 
the transition rates |12) -o- and |8) -f-^ |9) equalise, and we 
may transfer the coherences between the former to the latter, 
with a relative phase shift of n. We use wi = 200 MHz. 



(a) IWilO>p 




10 15 20 25 30 



(b) 



- - - \{m\'' 

|»|^(!|10>+|11»|- 




FIG. 4: Shows that at the u)q = 5 resonance, second order, 
two-photon transitions may be exploited since a; (10 — 9) ~ 
cj(ll — 10) [l^. A linear oscillating microwave field of strength 
uji = 200 MHz is used, (a) shows that driving at resonance 
prepares |10) —5- -i=(|ll) — |9)). The process is very sensi- 
tive to detuning from resonance, (b) and (c) illustrate how 
slight detuning of the microwave frequency may be used to pre- 
pare other superpositions such as |9) — )■ (j|10) + |11)) and 

|11>-^T75H|10) + |9». 



leaving only the e_5 term. Here i?_5 « — ^, so its eigen- 
value lies almost exactly half-way between the |±, —4) state 
eigenvalues: states |9), |11) have energies E± w i?_5±i?_4 
[ISj . This gives the striking feature at 2.3 GHz in Figl^I^b) 
where the |10) |9) and |11) o |10) lines coincide and 
where an EPR pulse would simultaneously generate coher- 



Hq = ilsS, + ujjlz + ASz (g) + BSz ® Ix- 



(18) 



Here £7^ = — w is the detuning from the external driv- 
ing field and w/ = Sluq is the nuclear Zeeman frequency. 
At resonance, Q,s = 0. As the hyperfine terms are weak, 
terms like Sx®Ix + Sy® ly are averaged out by the rapidly 
oscillating (microwave) driving. The remaining Hamilto- 
nian uilz -t- ASz ® Iz + BSz ® Ix conserves m^. For a spin 
S ~ ^,1 = ^ system like Si:P, the Hamiltonian decou- 
ples into two separate 2x2 Hamiltonians H^^^j.i. In the 
nis = +^ subspace. 



H. 



1 , A. 



B 



(19) 



where the Pauli matrices are defined relative to the basis 
\ms)(E)\mi) = \ + ^)®\zL^) while in the = — i subspace: 



1/ ^^ B 



(20) 



where the Pauli matrices are defined relative to the basis 
|tos) (g) |m/) = I - i) (g) I ± i). It is easy to see from Eq. ipO)) 
that if 



(21) 



only the ^ctx term remains. This is the "exact cancel- 
lation" condition. While reminiscent of hyperfine cancel- 
lations resonances, there are key differences, other than 
the obvious one that the former occur as a series when 
m = —uiQ in Eq.(IS]), while the "exact cancellation" occurs 
at a single magnetic field. In particular, since the EPR 
cancellation resonances of Eq.(l5]) affect both nuclear and 
electron spins, at m = — wqi the eigenstates assume "Bell- 
like" form: 



l*±) = 



1 

3 



(22) 

where the e, n subscripts have been added for clarity, to 
indicate the electronic and nuclear states respectively. In 
contrast, for exact cancellation, they give superpositions of 
nuclear spin states only: 



I*) 



1 



2>n ± I 2' 



(23) 
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which still permits interesting manipulations of the nuclear 
spin states [l6| . 

Note that, while exact cancellation eliminates the full 
Ising term ASz <E) Iz, the EPR cancellation resonance elim- 
inates only the non- isotropic part. In fact the EPR can- 
cellation resonances result in a variety of different types of 
resonances. As discussed below, the ujq — 7 A resonance 
does not cancel the hyperfine coupling at all; it equalizes 
the Bloch vector of the states in adjacent m-subspaces pro- 
ducing another effect (the "cancellation" label is in this 
case applied in a somewhat different sense). 

EPR cancellation resonances are in practice a much 
stronger effect than exact cancellation: decohering and per- 
turbing effects of interest in quantum information predom- 
inantly affect the electronic spins, not the nuclear spins. 
Exact cancellation appears in the rotating frame Hamilto- 
nian (which contains only terms of order MHz). It will not 
survive perturbations approaching the GHz energy scale. 
The cancellation resonances, on the other hand, arise in 
the full Hamiltonian, eliminate large electronic terms and 
can potentially thus reduce the system's sensitivity to ma- 
jor sources of broadening and decoherence. 

It is valuable to recall a major reason why the "exact 
cancellation" regime is so widely exploited in spectroscopic 
studies. In systems with anisotropic coupling, the spectra 
depend on the relative orientation of the coupling tensor 
and external field. Thus for powder spectra, which neces- 
sarily average over many orientations, very broad spectral 
features result. At exact cancellation, the simplification of 
the Hamiltonian is dramatically signalled by ultra-narrow 
spectral lines [l^ . 

III. SI:BI AS A TWO-QUBIT SYSTEM 
A. Basis states 



Adiabatic state 


high-field state 


logical qubit 


|12> 


l + i,-3i) 




n 1 \ 


n 1 \ 


l + i-4i) 




|leO„) 


q^i 

K/ 


1 - i -si's 




IOrl„) 


|iu; 


2' ^2/ 




\OeOn) 


Oontrolled operation 


Si:P transitions 


Si:Bi transitions 


Rv.ie), ® |o){o|„ 






<^ll-10 


Rv_{e)e ® |l)(l|n 


L04-1 




1^12-9 


\0){0\e(SRv_{0)n 


UJ2-1 




1^10-9 


\i}{i\e(»Rv.{e)„ 


(^4-3 




<^12-ll 



TABLE I: conditional single-qubit rotations of angle 6 about 
vector v_ in the Bloch sphere, denoted RJJ)), and correspond- 
ing transition frequencies. Frequencies in boldface correspond 
to qubit operations which are EPR-allowed at B = 0.1 — 0.6 
T; i.e. they require only fast (ns) EPR pulses. AH four EPR 
operations are possible for Si:Bi, whereas for Si:P nuclear qubit 
operations require slow (^s) NMR pulses. This scheme allows 
for cheap, controlled qubit operations, whereas single qubit op- 
erations would require twice the number of pulses. 

NMR transitions. Measurement of the qubits in the com- 
putational basis has to be performed at high fields, where 
the adiabatic logical qubit coincides with the electron and 
nuclear spin states. All simultaneous nuclear and elec- 
tronic qubit flips are forbidden for systems with isotropic 
hyperfine coupling A including both Si:P and Si:Bi. We 
note that, in spin-systems with "exact cancellation" and 
anisotropic A, the AI^ Sz coupling does permit simul- 
taneous nuclear-electronic qubit flips. These were recently 
shown for the organic molecule malonic acid [l6j ; the disad- 
vantage here is that single nuclear qubit rotations (essential 
for quantum computation) are not EPR-allowed. 



The adiabatic eigenstates of the well-studied 4-state 5* = 
1/2, J = 1/2 Si:P system can be mapped onto a two-qubit 
computational basis: 



Adiabatic state 



logical qubit 




With a 20-eigenstate state-space, the Si:Bi spectrum is 
considerably more complex. However, we can identify a 
natural subset of 4 states (states |9), |10), |11) and |12)), 
which represents an effective coupled two-qubit analogue. 
As hyperpolarization initialises the spins in state |10) 
and this state has both the electron and nuclear spins fully 
anti-aligned with the magnetic field, although it is not the 
ground state, it can be identified with the |OeO„) state. The 
other states - just as in the Si:P case - are related to it by 
adding a single quantum of spin to one or both qubits: 

For both systems, there are in principle 4 transitions 
which would account for all possible individual qubit oper- 
ations, as listed in Table H] We show below that for Si:Bi, 
all single qubit operations are EPR-allowed for B < 0.6 T. 
For Si:P, only the electronic qubit-fiips (the first two) are 
EPR-allowed; nuclear rotations require much slower, ^s, 



B. Universal set of quantum gates 

It is known that for universal quantum computation it 
suffices to be able to perform arbitrary single qubit rota- 
tions and a two-qubit gate such as the CNOT [131 . We now 
show how we may exploit the strong hyperfine interaction 
of the Si:Bi system to achieve this using only EPR pulses. 

Control of the electron spins is facilitated by the Hamil- 
tonian H = Hq + Vx/y(t) where V^/yit) = uJi cos{Lut)Sx/y 
represents the external microwave frequency field oscillat- 
ing along the x or y-axis. This may be written as: 



{cos{Lut)Sx/y — sin{u!t)Sy/x) (24) 



We label the first component the right handed (RH), and 
the second term the left handed (LH) rotating fields. In the 
rotating frame between two eigenstates |e) and \g) which 
satisfy the selection rule {e\Sx/y\g) = jy > 0, the Liouville- 
von Neumann equation for the reduced two-level system 
is: 



(25) 
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a ioi << Qq where fio = l-^e — Ag| is the energy difference 
between the two eigenstates. For most possible transitions, 
r] is positive and the resonance condition is satisfied by the 
RH component of the oscillating microwave field. But for 
transitions | — ^ | — ,m— 1), rj = a^b^_^ < and 
as such it is brought to resonance by the LH component 
of the microwave field. This feature, which is explained 
in more detail in AppcndixEl may aid in selective qubit 
manipulation as will be explained below. 

The transition strengths given in Eas. ([T^ . (jl3p and 
are given by \r]\'^. Eg. ((25)) shows that the qubit rotation 
speed, given a fixed microwave field strength, is determined 
by the mixing factor rj. As B —> oo, rj 1 for high-field 
EPR transitions, and ?y for high-field NMR transitions. 
At magnetic fields where A ^ B, however, mixing occurs 
and 77 will be appreciable for said transitions. 

At the TO = —4 cancellation resonance, established when 
ujQ ~ 4: A (B = 0.21 T), 1 77 1 for both nuclear and electronic 
qubit operations become exactly equal: this is simple to 
verify from Eas. (|T^ - ([ni) by setting 6^5 = and 0_4 = 
tt/2. We show numerically in Figl5lja)-(b) that this means 
a TT pulse on the nuclear qubit becomes as short as on the 
electronic one. For higher fields, however, 77 for nuclear 
qubit flips falls off as 1/B. 

The EPR pulses at our disposal allow us to perform con- 
trolled single qubit unitaries Rv{9) where v lies in the x — y 
plane. Two orthogonal Paulis suffice to generate arbitrary 
single-qubit unitaries [isj using at most three pulses, and 
we may construct the controlled and Hadamard gates 
by these pulse sequences: 



5/4. This ensures that at time 



(26) 



The possible controlled operations are shown in Table H] 
Single-qubit gates would require us to repeat the set of 
controlled EPR pulses for both the controlling qubit basis 
states, and as such would require twice the time. 



C. Selective qubit gates 

Consider the initial state IV') — (|lclri) + |Oeln)). If 
we wanted to perform a CNOT gate on this state, with 
the electron qubit as the control, we may choose to use the 
transition frequency a;i2-ii as dictated by Table HI How- 
ever, this frequency is only a few MHz different from that 
of wg-s, and a short pulse of 50 ns would have a band- 
width of ^ 20 MHz which would also drive the transition 
between states |9) and |8), and thereby effect a noise oper- 
ation on our qubits. There are two strategies to overcome 
this complication: 

1. by tuning the microwave frequency to be exactly be- 
tween the wanted and unwanted transition frequen- 
cies, we ensure they both suffer the same effective 
uj[ ^ wisincd + n^^^)), where T is the 
pulse duration and 51 is the transition frequency, 
such that the only variable affecting the two tran- 
sition rates would be the mixing factor ?/. Near 
the TO = — 4 cancellation resonance at B = 0.22 T 



we get 1^12^ 
t ~ , I — 1, we perform the operation: 



Ha^|lel„) + l|Oel„)) - (-*|leO„) + |0el„)) 



(27) 



This is shown numerically in FigjSJc) 



2. The transition |12) |11) utilises the RH component 
of the microwave field, whereas the |9) — ?> |8) transi- 
tion uses the LH one. By generating a RH circularly 
polarised microwave field, we would be able to select 
for the desired transition, as shown in FiglSfd). 



(=)•»„-,. 




B = 0.21 T 


---K*|0,o„)p 
— IW|i=o„)p 
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FIG. 5: (a) and (b) show Rabi oscillations utilising a linearly 
oscillating microwave field of strength L0i/2n = 200 MHz. For 
B — 0.21 T, the time taken for the electron qubit flip |OeOn) — >■ 
|le0n) is identical to that for the nuclear qubit flip |0e0n) — >■ 
|0eln). (c) and (d) show selective nuclear qubit flips, and have 
microwave fields of strength uji/2n = 100 MHz. This is because 
the energy difference between the eigenstates is smaller than 
in the previous two cases and lji must remain perturbative. 
(c) Utilises a linearly oscillating microwave field, which is non- 
selective for short pulses. At B = 0.22 T the rotation speed ratio 
is 7-^^"^^! = |. A Stt rotation of the nuclear qubit corresponds 



to a 47r rotation of the unwanted |9) — > |8) transition. In (d) 
we use a RH circularly polarised microwave field, which selects 
for the desired conditional nuclear qubit rotation, and is much 
more efficient than the linearly polarised case. 



D. Scaling with controllable Heisenberg interaction 

So far we have only described how to perform two-qubit 
gates in a single site of Si:Bi (or any other nuclear-electronic 
system obeying the same Hamiltonian and with a large 
enough hyperfine exchange term). This is very limited, 
however, and we need to be able to scale the system so as 
to incorporate arbitrarily large numbers of qubits. With 
Si:Bi, the possibility does exist to further utilise the 20 di- 
mensional Hilbert space, which would provide a maximum 
of 4 qubits. This is, however, not scalable, and we will 
still be limited to just 4 qubits as we cannot create more 
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energy levels within the single-site Si:Bi system. The only 
feasible option that remains, is to have spatially separated 
Si:Bi centers, between which we can establish an interac- 
tion. The original Kane proposal [l| envisaged a Heisen- 
berg interaction between nearest neighbour electrons which 
could be controlled via electric fields. This interaction has 
the form: 



H,nt = JS'.S'+i (28) 
As shown by [Tgj such an interaction can be used to produce 



a VSWAP gate: 



to spin-bath dynamics. Temperatures would also be main- 
tained at low levels so as to minimise the phonon-bath 
induced decoherence. Here, we wish to employ a phe- 
nomenological model of decoherence for nuclear-electronic 
systems, resulting only from stochastic magnetic-field fluc- 
tuations. Taking the Hamiltonian from Eq.(IT]) and adding 
to it a pcrturbative term involving independent temporal 
magnetic-field fluctuations in all three spatial dimensions 
(all of which take a Gaussian distribution with mean and 
variance a^) gives, in the interaction picture: 



H{t) = J2^nit)Sn{t) + UJ,,{t)Sinit) 



(33) 



VSWAP = e''te- 



(29) 



Using two such gates, together with single-qubit unitaries, 
we can establish a CZ gate between the electrons: 



VSWAPe"**- 1 VSWAP 
(30) 

As stated in previous sections, we cannot perform rotations 
about the z axis of the Bloch sphere directly, but we can 
use our EPR pulses about the x and y axes to produce the 
required single-qubit unitaries. 



e'^'e-'^y^e-'^'- Te- 



rn 



The CZ gate can be turned into a CNOT gate by sim- 
ply applying a Hadamard on the target qubit before and 
after the application of the CZ. Such an interaction affects 
the electron spin basis states, and not the adiabatic ba- 
sis states to which we have designated our logical qubits. 
Therefore, we must apply our electronic two-qubit gates in 
the high-field limit where mixing is suppressed, and where 
there is a high fidelity between the adiabatic basis and 
spin basis. A consequence of this is that the energy differ- 
ence between different eigenstates will be very large and, 
as is well known if \Ei — Ei^i\ >> J, the Heisenbcrg 
interaction between the two effectively becomes an Ising 
interaction JSl ^ S^. To use the above scheme of produc- 
ing entangling two-qubit gates between all 4 eigenstates in 
each of the two adjacent sites, we would have to establish 
a very strong J. 

Alternatively, we can set Bi and S^+i to be sufficiently 
different, and J small enough, such that we only get an 
Ising interaction between all relevant eigenstates. It is in 
fact easier to produce the CZ and CNOT gates with an 
Ising interaction, as it only requires one exchange operation 
and not two as in the case of the Heisenberg interaction 

MM. 



cz^ 



(32) 



IV. DECOHERENCE FROM TEMPORAL 
MAGNETIC FIELD FLUCTUATIONS 

For practical quantum information processing in silicon, 
the substance would need to be purified so as not to con- 
tain any ^^Si such that no decoherence would result due 



We may ignore the nuclear term as its effects will be neg- 
ligible. We then follow the standard procedure of deriving 
a Born-Markov master equation [2^ : 



dt 



{p{t))^^\{p{t)),H^+alY,Y. 



e X 



simp{t))s^{n) - - [{p{t)),si{n)s,,{n)\^]{M) 



where Sn{^) are the electron spin operators in the eigen- 
basis of Hq^ f2 is the energy difference between two such 
eigenstates, and /„ = (^^)~^- (p(0) is defined as the 
density operator of the coupled nuclear-electronic system, 
averaged over either an ensemble of such systems, or re- 
peated experiments on a single system. Further details for 
the derivation of Eq. (p4)) can be found in Appendix |B] 

The rate term e^^"^^ imposes the results of the adia- 
batic theorem into our master equation. The quantitative 
condition for adiabatic evolution is often cited as |2a| 



mum 



« 1. 



For the model described here, this translates to 



(35) 



{<P\Sn\i^) 



dB„ 
dt 



{E^ - E^,f 



{'P\Sn\i^) 



1 



« 1 (36) 



which means that, in the case of 



\Sn\^P) > 0, if the 
magnetic field is fluctuating sufficiently slowly, the proba- 
bility of transition between eigenstates ji/') and |0) becomes 
vanishingly small. 

Since we have made the rotating wave (or secular) ap- 
proximation, and are only interested in the interaction pic- 
ture dynamics of our system, we may drop the Hamilto- 
nian commutator in Ea. (j34p . leaving only the dissipator 
term. We are now equipped with the tools to address 
decoherence in our quantum system. We may model our 
gates as ideal unitaries that can prepare some superposi- 
tion IV') = a\e) + b\g) between the adiabatic basis states 
|e) and \g), which is then decohered according to our noise 
model. The dephasing and depolarising rates are deter- 
mined by applying our master equation and measuring the 
rate that the observables 



^Jtr{a,p{t))^ +tr{ayp{t))^ 



(37) 



9 



and 



(38) 



decay respectively. In the cases that these decay exponen- 
tially, we may characterise the dephasing and depolarising 
times by T2 and Ti respectively, which are the inverse of 
the decay rates. Such times are measured in EPR exper- 
iments [2^, 111]- The Pauli matrices denoted here are in 
the eigcn-basis of the reduced two-level system in question. 
We will study two types of noise; Z noise and X noise, so 
named due to Gaussian magnetic-field fluctuations in the 
z-axis and x-axis respectively. Throughout this section, 
when an adiabatic state is designated with m, it is implied 
that \m\ < {I + i), and states where m = ±(/ + i) are 
explicitly designated. 



Z noise 



(a) p(to) 




,„ (c) X noise B = 6 T 

.p(io+2.0/a2) 



20 



20 

(d),Z noise B.= 0.188 T 
5 0.5,-- ■ ,i(io +20/q2) . 
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20 

■ ,(e) Z noise B = 6 T 
>(.io+20/a2) 




20 



20 



FIG. 6: Shows the change in density operator elements of — 
-^|9) + -^112) after a period of , given diabatic Gaussian 

noise of variance a^, at two field regimes: the frequency minima 
of B = 0.188 T, and B = 6 T . (a) Population of p in the 
eigen-basis of Ho at time to- (b) - (c) show the effect of B on 
the properties of X noise, (b) Near the frequency minima, X 
noise couples every part of the full Hilbert space, and effects a 
depolarising channel for the full system, ultimately resulting in 
2^1. (c) At B = 6 T, X noise decouples states |-|-, m), | — , m — 1) 
from the rest of the Hilbert space, and effects a depolarising 
channel in that subspace only, (d) - (e) show the effect of B on 
the properties of Z noise. Z noise conserves angular momentum 
and hence keeps to the 4-dimensionaI Hilbert space of m = 
—3, m — 1 = —4. (d) shows that at the frequency minima of 
B — 0.188 T, Z noise effects independent depolarising channels 
for each m subspace. As a result, the population of states |12) 
and |9) equalise with those of states |8) and respectively. (e) 
shows that at B = 6 T, {+,m\Sz\ — ,m) ~ and we simply get 
a dephasing channel for the | + , ?n), | — , m — 1) subspace. 

Given Z noise, we may consider our system as a decou- 
pled four-level system with sub Hamiltonian: 



Hsub = {^mCTx + AniCr^) (flm-lCrx + '^m-lCTz) (39) 



ignoring any identity terms. This is possible as there will 
be no transfer of population to other components of the 
Hilbert space. We may write Sz in the adiabatic basis 

{|+,m), |-,to), |+,m- 1), |-,TO- 1)}: 



/ cos(6',„) -sin(6'„i) 

-sin(6'„) -cos(6',„) 

V 






cos(6'„_i 
- sin(6'„,_: 



Substituting this into Ea. (|M)) gi 






-sin(6'„_i) 
) -cos{9,n-i)J 



(40) 



n— m— 1 
^,2 



;ives 



nr. 



^cos2(0„)(a;'(p(t))ar-(p(i))) + 



{\+,n){-,n\{p{t))\-,n){+,n\ - 
\+,n){-,n\\^,n){+,n\{p{t)) + H.C.).{41) 

In the infinite-field limit our EPR local unitaries can only 
create superpositions a\+,m) + b\ — ,m — 1). As the noise 
operator takes the form (T™ © this regime, this 

superposition may be considered to exist as a decoupled 
two-level system. We may therefore solve Ea. (|4T|) (for the 
two-level subspace in question) analytically: 



(42) 

where C is the Liouville superoperator, whose action on p 
is given by Eq. (PT|) . This is simply the dephasing channel 
for a spin 1/2 particle [l8| : 



£{t) o p = {1 - X(t))p + X{t)azpaz 



(43) 



with probability X{t) of performing a operation under 

_ t 

conjugation. Here, X{t) — ■^—^ with a T2 time of 

This is illustrated by FigEl^e), where only the off-diagonal 
elements of p{to), as shown in FigElJa), decay. At low fields 
however, m 15*2 1 — , m) > and we cannot ignore the 
exchange term in Ea.(pT|). In this case, we may use the 
4-dimcnsional Bloch vector representation of our density 
operator: 



if ' \ 

p{t) = ^ E "^^3(^)^1 «) (Tj , noo(t) = 1 

\ij=0 / 



(44) 



It is possible to map the dynamics of the density operator 
to that of the Bloch vector [23 as 



dn{t) 
dt 



Cn{t). 



(45) 



For the stated problem the 16 simultaneous differential 
equations can be solved to obtain analytic expressions for 
the dephasing and depolarising rates. Alternatively, by 
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decomposing n{t) in the eigenbasis of £, denoted n.; with 
generaUy complex eigenvalues A/ , we may represent the dy- 
namics of the Bloch vector as 



1=0 



(46) 



where ci are determined by the initial conditions. It is 
the real component of the eigenvalues which leads to decay 
in population of the eigenstate. The infinite time state 
is therefore a superposition of eigenstates n.; such that 
Re{Xi) = 0. 



1. Adiabatic Z noise 

Here we may set / — oo =4> e~^^^ ~ for VP' > 
0. There will be no dcpolarisation in this case, and we 
may only have pure dcphasing. For superpositions of type 
a\+,m) + 6| — ,TO — 1), the dcphasing rate, paramcteriscd 
as the decay of the off-diagonal elements of the subspace 
in question, is given by 



1 a ( I 



(47) 



When cos(0„i) = — cos(0m-i), which is satisfied at the 
frequency minima, dcphasing due to Sz is completely re- 
moved. 

For superpositions of type a|±,m) H- 6|±,to — 1), the 
dcphasing rate is given by 



1 



7f- = -^icos{0rn) - C0s(6'„_i))^ 



(48) 



Here there are two regions where the 5"^ caused dcphasing 
is removed; when cos(6'm) = cos(6',„_i) which occurs at 
the frequency maxima, and at the high-field limit where 
cos(0„i) — cos(0„i_i) — 1 \/ m, rendering such transitions 
as only NMR allowed. 

For superpositions of type a|±, ±(/ + i)) + b\zL, m) the 
dcphasing rate is given by 



1 

7^ 



(49) 



2. Diabatic Z noise 



Here we may set / — )- 



ing the Bloch vector differential equations yields analytic 
expressions for the dcphasing rates. For an initial super- 
position of a\+, m) -\- 6| — , m — 1) this gives: 



1 



(cos(6'„,)cos(6'„,._i) + 1) 



and for a|±, m) 4- 6|±, m — 1) : 



1 

T'2 



(cos(6'„)cos(6',„_i) - 1) 



(52) 



(53) 



Eq. (|52p reaches a minimum value (hence giving the longest 
T2 time) when cos(0„i) = — cos(6'm-i) i-e. at the frequency 
minima. Unlike the adiabatic Z noise case, this value does 
not reach 0, but rather reaches approximately half its max- 
imal value at the high- field regime. Conversely, Eg. ((53)) 
reaches its maximal value when cos(0m) = — cos(0m_i), 
attaining approximately the same value as for Ea. (j52p at 
this regime. Note that unlike the adiabatic case, there is no 
decoherence minimum at the frequency maxima; the decay 
rate simply vanishes as S — >■ 00. 

Because the exchange terms in Eq. (|lT|) contribute to the 
dynamics for diabatic Z noise, there will also be depolaris- 
ing noise in each m subspace, equalising the population in 
states |±,m) . FiglHlJd) shows the effect of this dcpolarisa- 
tion at the frequency minima. This dcpolarisation rate of 
each m subspace is given by: 



1 



■ sm 



(54) 



which vanishes as i? — ^ cxj, and maximises at the avoided 
crossing cancelation resonance. Given any superposition 
\fPg\9) + e*'^A/^|e), with states \g) and |e) each existing 
in a different m subspace, such that > Eg, the dcpolar- 
isation is given by: 



-P. I 



-tjTl 



-tin 



) (55) 



where 1 jT^^'^ is the dcpolarisation rate in the m subspace 
for |f;) and |e), respectively. The steady state solution for 
diabatic Z noise given such superpositions is given by 



which reaches its minimal value of as B — > cxj, whereas 
for superpositions of type a|±,±(/ -I- i)) -I- 6|=F,™) it is 
given by 



1 a' 

^ = T^°^ 2 



(50) 



which reaches its minimal value (which is generally greater 
than 0) at i? = T. The steady state solution for adiabatic 
Z noise is given by 



n(oo) = 1 (K) l-f cil (g) (T^ + C20z «) 1 + csCT^ ® Oz- (51) 



n(oo) = 1 (g) 1 ci(T2 ® 1 



(56) 



where c\ G [1,-1]. 

For superpositions of type a|±, ±(/ + i)) + 6|±, m) the 
dcphasing rate is given by 



(1 - cos(&^)) 



(57) 



and for superpositions of type a|±, ±(/+ i)) + fe|=F, m) the 
dcphasing rate is given by 



1 , 

^ = Y(l + cos(e,n)). 



(58) 
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The depolarisation rate can be calculated as in the previous 
case, using Eq. ([55|) . and noting that one of T^^^ is equal to 
oo. The steady state solution for diabatic Z noise for such 
superpositions is given by 



n{oo) = 1 1 + cio-^ (g) 1 + C2I ® (7^ 



(59) 



For diabatic Z noise, Fig (7] shows the analytical depolarisa- 
tion and dephasing rates for subspace m = — 3,m— 1 = —4. 
FiglHa) shows the numerically calculated T2 times for all 
EPR lines depicted in FiglUJb), whilst FiglHl^a) shows the 
numerically calculated depolarising times Tl in units of 
for each m-subspace. 




a.| + ,-3> +fj|-,-4) 

---o|±,-3> + 6|±,-4) 



(b) 



B / T 



B / T 



FIG. 7: The exponential decay rate given by F in units of 
/2 for diabatic Z noise driven (a) depolarisation and (b) 
depliasing. This is done in the 4 dimensional subspace of 
m = —3, m — 1 = —4. (a) shows that in each subspace, the 
depolarisation rate maximises when 6m = 7r/2, or the avoided 
crossing cancelation resonances, (b) shows that at the high-field 
limit, the dephasing rate of a\-\-,m) + h\ — ,m — 1) is maximal, 
whilst that of a|±, m) 4- 6|±, m — 1) becomes vanishingly small. 
These rates both approximately reach the value of 1/2 at the 
frequency minima. 



B. X noise 



(a) Z noise 
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(b) X noise 
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FIG. 8: Simulated dephasing times in units of 2/o? for dia- 
batic (a) Z noise and (b) X noise. In (a) the superpositions 
a\+,m) + b\-,m- 1) have T2 times of at B > 0.6 T, and 
approximately 4/a^ at the frequency minima. Superpositions 
a|±,m) + 6|±,m — 1) also have T2 times of 4/a^ at the fre- 
quency minima. However, as B increases, these become NMR 
transitions and will have T2 times of 2/{cP5). The colour bar 
has been truncated after 3 to aid visibility. In (b), the T2 time 
does not vary by much, and reaches its maximal points at fields 
less than the frequency minima. 



(a) Z noise 



(b) X noise 
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X noise is less trivial, as it couples all components of the 
Hilbert space so we cannot consider a sub-Hamiltonian in 
isolation. Solving the resulting 400 Bloch equations (for 
Si:Bi) would be unfeasible, so only numerical calculations 
are given here. Furthermore, the adiabatic condition must 
be violated for X noise to have any effect, as there are no 
S^{VL = 0) terms in Eq.(l34l). In the high- field hmit the 
X noise operator will take the form of ®t in the basis 
{|-f , m), | — , TO — 1)}. At such fields, as shown in FiglHlJc), 
an arbitrary superposition of a|-|-, to) — , to — 1) suffers a 
two- level system depolarising channel. At low fields, how- 
ever, the dissipation is not contained within the to, to — 1 
subspace, and as indicated by FiglHlJb) the system eventu- 
ally decays to ^1. 

For X noise, all dephasing is a result of the depolarising 
noise that is effected by the X noise operator, and as shown 
in FiglSJb), at B > 0.6 T the dephasing time is A/a^ for aU 
transitions. This value increases only slightly at magnetic 
fields smaller than the frequency minima for transitions 
involving m < 0. 

FiglHJb) shows the different forms of depolarising rates 
for X noise. At high magnetic fields, where the only non- 
vanishing matrix elements of the X noise operator are 



FIG. 9: Simulated depolarising times for diabatic Z and X noise 
with Qf^ = 9 MHz. (a) Given Z noise, the decay of tr{a,p{t)) 
within each m-subspace is always exponential. For m < the 
Tl time reaches a minimum at the avoided crossing cancellation 
resonances. The Ti time for subspace m and m — 1 become 
identical at the frequency minima (b) Given X noise, the decay 
in the | — ,m — 1) subspace follows an exponential curve 

at high magnetic fields, but at low magnetic fields such as the 
frequency minima, it follows a double exponential fit. 

m\Sx\~,m— 1), the depolarising noise follows an expo- 
nential decay. Under intermediate magnetic fields however, 
the dissipation follows a more complicated mechanism and 
the decay is better explained by a double exponential fit. 

V. CONCLUSIONS 

A coupled nuclear-electronic spin system with large A 
and / > 1 will give cancelation resonances at non van- 
ishing magnetic fields. In the intermediate-field regime 
{B ~ 0.1 — 0.6 T) the exceptionally large values of A and 
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/ of Si:Bi generate a series of such cancellation resonances. 
Two eigenstates, with an EPR transition, which have a 
cancelation resonance will have two ^ = points in the 
spectrum; a frequency minimum and maximum. These 
points arc associated with subtle and hereto unstudied ef- 
fects: both line broadening and decoherence effects may 
be reduced. If the electronic and nuclear spins of Si:Bi are 
used as a coupled two-qubit system, the cancellation res- 
onances allow a universal set of quantum gates to be per- 
formed with fast EPR microwave pulses, eliminating the 
need for slower radio frequency addressing of the nuclear 
qubit. 

One scheme would envisage the following stages: (1) hy- 
perpolarization of the sample into state |10) (in which the 
2-qubit system is initialized as |0)„|0)e) at i? w 5T (2) A 
magnetic field pulse (~ 10 T/ms, of duration lower than 
decoherence times) would reduce B to ~ O.IT (3) A se- 
ries of EPR pulses would be carried out on the system 
in the low-field regime (4) As the magnetic pulse decays, 
the system is restored to the high-S limit, leaving it in 
the desired superposition of |0)e|0)„, |l)e|l)n, |0)e|l)n and 
|l)e|0)n basis states. (5) Turn on a Heisenberg exchange 
interaction between spatially separated bismuth sites to ef- 
fect a two-qubit gate on the electron qubits. Thus, given 
the capability to rapidly (< 1 ms) switch from the high 
to intermediate-field regime, Si:Bi confers significant addi- 
tional possibilities for quantum information processing rel- 
ative to Si:P. We conclude that in particular the m = —4 
cancellation resonance at 0.21 T may represent an "opti- 
mal working point" for quantum computation with Si:Bi 
for several reasons: (a) it is the point at which nuclear 
and electronic qubit flip times equalize (b) it provides 
df/dB = points (c) it is associated with line narrowing. 
The other resonances provide a subset of these character- 
istics. The advantages of the df/dB = points have been 
quantified relative to Gaussian noise: the Markovian model 
of decoherence presented here suggests longer dephasing 
times at the frequency minima given Gaussian temporal 
magnetic field fiuctuations in the z-axis. 
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Appendix A: Selective rotations 

Consider the coupled nuclear-electronic spin system in 
the eigcnbasis of the Hamiltonian Hq given by Eq. ([T|) . The 
electron spin operators in this basis are given by the unitary 
transformation 



S' =V''S^V S', 



V^SyV 



(Al) 



where y is a matrix whose i*'' column is the i*^ eigenvector 
of Hq. We want to be able to isolate two eigenstates of 
this Hamiltonian, and perform unitary dynamics in that 
subspace. Tracing out all eigenvectors other than |e) and 

\g) gives 



■r \eg 



1 







2 \\g){g\iS+ + S^)\e){e\ 



\e){e\iS+ + S^)\g){g\ 




2 l^|.g)(.g|(^_-^+)|e)(e| 



\e){e\{S^-S+)\g}{g\ 




(A2) 



where is a measure of the degree of mixing of the spin 
basis states in the eigenstates of the Hamiltonian. As the 
absolute energies given by the eigenvalues are meaning- 
less physically, we can re-scale the eigenvalues of Ha by 
adding to it an identity term such that Ae and 

Xg are the eigenvalues of eigenvectors |e) and \g) respec- 
tively, where Ae > Ag. This gives: 



^0 ^ m = 



\g) 



rjrem 
^0 



(A3) 
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such that rio = I'^e ^ Ag|. Given a perturbative Hamilto- 
nian of the form in Eg. ip^ . and assuming that all EPR- 
allowed transition frequencies are unique, we may solve the 
Liouville-von Neumann equation for the two-level subsys- 
tem in the rotating frame of Hq: 



.^1 



p(t), e'^**^ { [cosinot)iS',/^r + sminot)iS'^/,r) + (cos(l]ot)(^;/,)^^ - sininot) {S'^/,r) } , 



. fin . . 



= «^ [pit), ^x/y] + [p(i)) cos(2rioi)o-^/y - sin(2fioi)o-y/x] 



(A4) 



r 



There are two possible regimes for the dynamics of this 
system. Those for which t] is positive, and those for which 
r] is negative. By inspection of the eigenvectors of Hq given 



by Eqs.®,®, it becomes clear that the latter occurs only 
for transitions of type \ — ,m) O | — ,m — 1). Table HIl shows 
the form that matrices [S'^/yY^ take for both regimes; the 
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only one which undergoes a sign change for 77 < is the 
{S'^Y^ matrix. This has the effect of inverting the x-axis of 
the Bloch sphere. In the case that > 0, the circular polar- 
isation needed to achieve resonance with the Hamiltonian 
is of the form co&{VlQt)Sj./y+sv[\{Vli)t)Sy/x. We call this the 



right-handed (RH) field. When 77 < 0, the circular polar- 
isation must be of the form cos.(VtQt) S ^ / y — sm{^lt)t)Sy/x, 
which we call the left-handed (LH) field. 





right-handed 


left-handed 




> 


< 




1 ( \e){e\S+\g){g\\ |,|^ 
' \\g){g\S-\e}{e\ ) ^ - 


1 ( \e){e\S-\g){g\\ |,|^ 
^ \\g}{g\S+\e){e\ ) 


(s'yr 


, ( -|e>(e|S+|3)(3|\ _ |,|^ 
' \\g){g\S.\e){e\ j ^ ^ 


, / \e){e\S.\g){g\\ _ |,|, 
' \-\g}{g\S+\e){e\ J ^ y 








rotating field 







TABLE II: Pauli operators for the truncated 2-level system under resonance in the right-handed and left-handed regimes. The 
bottom row indicates the polarisation that the Vx/y{t) term must have in order to effect a e^'^'^x/y operator in the rotating frame. 



Appendix B: Master equation derivation 

Taking the Hamiltonian from Eq.(IT]) and adding to 
it a perturbative term involving independent temporal 
magnetic-field fluctuations in all three spatial dimensions, 
all of which take a Gaussian distribution with mean and 
variance a^, gives in the interaction picture: 

3 

= ^c^„(t)5„(i)+w„(t)«„(t). (Bl) 

n=l 

We may ignore the nuclear term as its effects will be negli- 
gible. We then write the Liouville-von Neumann equation 
in differential-integral form and take the average over the 
field fluctuations. Noting that {f^pit)) = ftiPi'^)) 

we may 

write this as: 



be another Gaussian function of the form 

{u;n{t)ujn{t - r)) = al-^e~^ (B3) 

where /„ = {^^) ^ is an inverse function of how fast the 
magnetic fleld fluctuates. In the limiting case of — > 
00, the correlation function tends to a'^S{T). To aid our 
calculation, we may write the noise operators in the basis 
of the eigenstates of Hq labeled |a),|6): 

n 

Sr^m = Y,^i^ba~ma){a\Sn\b){b\. (B4) 

a,b 

Moving such operators to the interaction picture simply 
gives e~'^*S'„(rJ). This gives: 



^ / dsK(t)L.„(s)) \\{p{t)),Si{t)] ,5„(s)j (B2) 



3 „t 
n=l •^*0 



where assigning Sn{t) = S'J(i) is valid as it is a Hermitian 
operator. {p{t)) is the density operator for the nuclear- 
electronic system, averaged either over an ensemble of such 
systems, or over many repeated experiments on the same 
system. Here, we have assumed that the field fluctuation 
statistics are independent of the quantum state of our sys- 
tem. These assumptions and approximations lead to a 
Born- Markov master equation. As U!n{t) follows a Gaus- 
sian distribution with mean 0, the flrst term of this equa- 
tion vanishes. Given that the correlation function drops to 
zero at flnite values, we may set the integration limits to 
to = and t ~ 00, and change the integration constant 
to r = t — s. As the temporal fluctuation takes a Gaus- 
sian distribution, we may set our correlation functions to 



dt 



(Pit)) = 



ri=l O.Sl' 



dT- 



(B5) 



If << meaning that the dynamic time scale of 

our system is much shorter than that of the decoherence 
caused by the magnetic-field fluctuation, which is a reason- 
able assumption for systems of interest, we may make the 
rotating wave approximation (often also referred to as the 
secular approximation) and drop all terms where ^ 57'. 
Furthermore, noting that [A, [B,C]] = ABC-BCA + H.C. 
(Hermitian conjugate) if A, i?, C arc Hermitian operators 
leads to : 



dt 



i—n ri 

+H.C. (B6) 



14 



We may decompose the integrand to 



1 



OO 

dr- 

2 V nfn 



e-fce'"" =7„(r!) + ir„(r!) (B7) 



where: 



-{p(t))=t[{p{t)),Ho + HLs 



1 r 



5t(l])(p(t))5„(0) - - [{p{t)),SimSnm\^] (B9) 



=e e 



n -(BS) 



Plugging these expressions into Ea. (jB6p . and moving back 
into the Schrodinger picture, gives us our final master equa- 
tion: 



where 



Hls = Y.Y.^nmslin)Snin) (BIO) 



n O 



is the Lamb shift and changes the energy levels of the sys- 
tem. This is a negligible effect and hence can be ignored. 
We use [A, _B]^ := AB+BA as the anticomutator operator. 
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